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("^ ' We present the foundations of a new approach to the Casimir effect based on classical ray optics. 

^^ I We show that a very useful approximation to the Casimir force between arbitrarily shaped smooth 

Cn , conductors can be obtained from knowledge of the paths of light rays that originate at points between 

(-H ' these bodies and close on themselves. Although an approximation, the optical method is exact for 

r^ ' flat bodies, and is surprisingly accurate and versatile. In this paper we present a self-contained 

derivation of our approximation, discuss its range of validity and possible improvements, and work 

out three examples in detail. The results are in excellent agreement with recent precise numerical 

analysis for the experimentally interesting configuration of a sphere opposite an infinite plane. 
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(3 , I. INTRODUCTION 

^^ , Revolutionary new experimental techniques have made possible precise measurements of Casimir forces 1]. Casimir's 

f| ' original prediction for the force between grounded conducting plates due to modifications of the zero point energy of 
Qh. the electromagnetic field has already been verified to an accuracy of a few percent. Variations with the conductor 
geometry and the effects of finite conductivity and finite temperature will soon be measured as well. Progress has been 
G ' slower on the theoretical side. Despite years of effort, Casimir forces can only be calculated for the simplest geometries. 
Beyond Casimir's original study of parallel platesl^, we are only aware of useful calculations for a corrugated plateQ 
and for a sphere and a platej^. The former was obtained with functional integral techniques quite special to that 
geometry and the latter was obtained by computationally intensive numerical methods. Simple and experimentally 
interesting geometries like two spheres, a finite inclined plane opposite an infinite plane, and a pencil point and 
a plane, remain elusive. The Proximity Force Approximation 5] (PFA), which has been used for half a century to 
%2^ ' estimate the dependence of Casimir forces on geometry, was shown by Gies et al. ^ to deviate significantly from their 
5t ! precise numerical result for the sphere and plane. Thus at present neither exact results nor reliable approximations 
are available for generic geometries. It was in this context that we recently proposed a new approach to Casimir 
effects based on classical optics |y|. The basic idea is extremely simple: first the Casimir energy is recast as a trace of 
the Green's function; then the Green's function is replaced by the sum over contributions from optical paths labelled 
by the number of (specular) reflections from the conducting surfaces. The integral over the wave numbers of zero 
point fluctuations can be performed analytically, leaving 

Here tr{x) is the length of the closed geometric optics ray beginning and ending at the point x and reflecting r times 
from the surfaces. Ar(x) is the enlargement factor of classical optics [3, 13, also associated with the r-reflection path 
beginning and ending at x. T>r is the subset of the domain, T), between the plates in which r reflections can occur. 
The factor (—1)'^ implements a Dirichlet boundary condition on the plates; different boundary conditions require 
different factors. Both l!,r{x) and Ar(a;) are very easy to compute either analytically in simple cases, or numerically in 
general. Ar(a;), although well known in optics, may not be familiar in the context of Casimir effects. We will describe 
its properties in some detail. 
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Eq. (jl.ip turns out to be a powerful tool to compute Casimir effects for generic geometries, and to identify, interpret 
and dispose of, divergences. Eq. (|I.1|I is not exact. Instead it is an approximation which is valid when the natural 
scales of diffraction are large compared to the scales that measure the strength of the Casimir force. In practice 
this will typically be measured by the ratio of the separation between the conductors, a, to their curvature, R. 
Although approximate, the optical approach is surprisingly accurate, as well as physically transparent and versatile. 
It generalizes naturally to the study of Casimir thermodynamics, to the study of energy, pressure, and momentum 
densities, to various boundary conditions, to fermions, and to compact and/or curved manifolds. This is the first 
in a series of papers intended to provide an introduction to the optical approach to Casimir physics. Here we will 
focus on fundamentals: how to derive the optical approximation and how to apply it to practical calculations of 
Casimir forces. In later papers in this series we study Casimir effects at finite temperature, the calculation of local 
observables like the energy density and pressure, and the generalization to conducting and other boundary conditions. 
Our first aim is to familiarize the reader with the use of the optical approximation, since this method of calculation is 
unfamiliar. In Section II we present some examples of the use of the optical approximation. First we review in more 
detail the treatment of parallel plates already presented in Ref. [g ■ Although it is no great triumph to rederive this 
classic result, the optical derivation illustrates several characteristic features of the method: rapid convergence, simple 
disposal of divergences and ease of computation, in particular. Next we present the case of a sphere and a plate. This 
too was summarized in Ref. 6|. Here we concentrate especially on the enlargement factor, both its interpretation 
and how to compute it. Also we illustrate the generic way that divergences can be eliminated. The numerical results 
we present here are more accurate than those of Ref. Ti]. Finally we apply the optical method to the case of a 
finite plate suspended above an infinite conducting plane - the "Casimir pendulum" . We show how all reflections 
can be computed and how the optical result differs from the proximity force approximation. In collaboration with 
O. Schroeder we are preparing a thorough study of the hyperboloid ("pencil point") near an infinite planej^]- In 
Section HI we discuss the derivation of the optical approximation from exact expressions for the Casimir energy. We 
show how a uniform approximation to the propagator turns into a uniform approximation for the Casimir energy. The 
derivation illustrates the nature of the approximation and shows the way toward improvements, which, in essence, 
amount to including the effects of diffraction. We present results for a massive scalar field in N dimensions in Section 
III. Higher spin fields will be considered in a later paper of this series. We discuss the general problem of divergences. 
The Casimir energy is generically divergent — or more properly, it depends in detail on the cutoffs that limit the 
conductivity of real materials at high frequency. However it is known that the Casimir force between rigid conductors 
is cutoff independent [1^. In the optical approximation the cutoff dependent terms in the Casimir energy can easily 
be isolated and shown to be independent of the separation between conductors. They therefore do not contribute to 
forces and can be dropped. Corrections to the optical approximation will bring in new surface divergences. In Section 
nil CI we discuss the relation of the optical approximation to previous works on "semiclassical" approximations to the 
Casimir energv[llj. In the last section we summarize our results, discuss their implications, and mention extensions 
to other interesting geometries. 

II. THREE EXAMPLES 

In this section we present three examples of the use of the optical approximation, cq. (|I.1|I . Our aim is expressly 
pedagogical: we want to demonstrate that this method can yield interesting and accurate results without onerous 
calculations. 

A. Parallel plates 

Casimir's original result for parallel plates can be derived in many ways. We present a derivation from the optical 
approximation in order to illustrate several generic features of the approach in the simplest possible context. The 
points we wish to stress are: ease of calculation; the rapid convergence in r, the number of reflections; and the simple 
and accurate treatment of divergences. The "semiclassical" method[ll| and the method of images'!^ generate exactly 
the same calculation as ours for parallel plates. However they do not generalize to less trivial geometries (although 
one might say that our method is the correct generalization of the method of images). We study a massless scalar field 
for simplicity, and quote the generalization to a massive scalar in a later section. For a flat surface the enlargement 
factor reduces to l/£r{x), so the contribution of the r reflection path is 

,,.-|,(-l).M„/^^A^, (ILl) 
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FIG. 1: Optical paths for parallel plates. The initial and final points on the paths, which coincide, have been separated so 
the paths can be seen, a) Even reflections 2, 4, and 6. Path 2' is distinct from 2 and illustrates the origin of M2„ = 2. b) 
Odd reflection paths. The paths shown form a family of continuously increasing length. Another family begins with the first 
reflection from the top. 



where Mr is the multipHcity of the path. It is convenient to separate the paths into "odd" (r = 2ri + 1) and "even" 
(r — 2n) according to the number of reflections. Some of these paths are shovirn in Fig. pp. Odd and even paths differ 
dramatically in their contribution to the Casimir effect: they differ in sign and in multiplicity: Mr = 1 for odd paths 
and Mr = 2 for even paths, as shown in the figure. The length of an even path depends only on n, whereas the length 
of an odd path varies with position. Finally, odd paths contribute a divergence to £, but do not contribute to the 
Casimir force. The even paths are finite and give the entire Casimir force. First consider the even paths. The length 
of the 2n refiection path is ^2n = 2na independent of x, as can easily be seen in Fig PP. The volume of each domain, 
X'2n, is the volume between the plates, Sa. Hence the contribution from even paths is 

£cvcn = --z-^^SaS^ — — TT = -7-777-^5'. (II. 2) 

27^2 ^-^ (2nar 14400^ 
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which is the famous result due to Casimir [J. Next consider the odd paths. There are two families. One is illustrated 
in Fig.^ The other family begins with the first reflection from the top plate. Their contributions are identical, giving 
an overall factor of two. The r — 2n + 1 reflection paths range in length from 2na to 2{n + l)a as can be seen from 
Fig. p], and contribute 

£2n+i^^2Sj dz^^, forn = 0,l,2,... (113) 

The first refiection contribution diverges at the lower limit. As discussed in the Introduction (and further in Section 
III) the divergence indicates dependence on the properties of the material composing the plates and is cutoff at a 
distance scale e determined by the microphysics. For example we can take e to be the skin depth or regard e as ~ c/A, 
where A is a frequency cutoff, for example the plasma frequency of the metal. Inserting e as the lower limit for n — 
and summing over n, we obtain the contribution of odd paths, 

he f°° 1 he 

£.,,^—,2SJ^ dzj^^^^-^^S (llA) 

This contribution displays the cubic surface divergence expected for a scalar field obeying a Dirichlet boundary 
condition [l3|. However, the divergent term — and indeed the sum of all odd refiections — is independent of a and 
therefore does not contribute to the force between the plates. Until now we have not considered the contributions 
from one-reflection paths that lie below the bottom plate or above of the top plate. It is easy to see that the sum of 
these contributions is identical to eq. (jll.4p and does not contribute to the force. This simple calculation illustrates 
some general features of the optical approach: 

• The even refiections dominate, give rise to attraction, and their sum converges rapidly in n. They are also 
attractive for Neumann boundary conditions, where the factor (—1)'^ is absent. They would be repulsive if 
one surface were Neumann and the other Dirichlet. In the case of parallel plates 92% of the Casimir effect 
comes from the second refiection, 98% from the second and fourth, and 99.3% from the second, fourth and sixth 
refiections. Similar results will be found to hold in more complicated geometries. 



• 



The only divergent contribution conies from the first reflection, ft does not depend on the separation and 
therefore does not contribute to the Casimir force. This result is quite general. To see the general argument, 
reconsider the first reflection from the bottom plate, Si 



£i,Sr = T^S / dz-—^ = —^S / dz-—^ - —^S I dzj^^-^. (II.5) 



The first term in eq. (jff.5|) combined with the contribution of the 1-reflection path outside of the plates (from the 
lower face of the bottom plate) is the cutoff dependent energy of an isolated plate. It is manifestly independent 
of the presence of any other conductor, and gives no contribution to Casimir forces. The second term is a 
finite effect of the first reflection. For parallel plates the finite contribution of the first reflection is cancelled by 
higher odd reflections. This occurs whenever the enlargement factor is l/l^, that is, when all the conductors 
arc planar. For non-planar surfaces the first reflection gives a (relatively small) cutoff independent contribution 
to the force. 

• The optical approach gives the exact answer for infinite plates. However it will fail when S^^"^ w a for the same 
reason that the capacitance of two finite, parallel metallic plates contains corrections of order a^/S'[lJ|: It is a 
poor approximation to consider the electric field inside two far separated plates (a ^ S"^") as constant inside 
and zero outside. Likewise, in the same limit it is a poor approximation to expect the Green's function for the 
field (j) to have contributions only from optical paths. The corrections, or edge effects, can be regarded as due 
to diffractive rays coming from the edges of the plates [T^ ■ We discuss corrections to the optical approximation 
in further detail in Section III. 

• The difference between even and odd paths has a fundamental origin, as already noticed in work on the "semi- 
classical" approximation to the Casimir energy jlll |. The even paths are truly periodic, in the sense that the 
momentum of the particle, aft er g oing around the path, returns to its initial value. These are therefore the paths 
that according to Gutzwiller 16\ contribute most to the oscillations of the density of states. The connection 
between these paths, the oscillation of the density of states, and the finite part of the Casimir energy has been 
noted many times[l7l| and is exact for parallel plates and related geometries {eg. flat manifolds with various 
topologies). However, the exactness of this result is an accident due to the particularly simple geometry. For 
example, there are very simple geometries in which periodic paths do not exist at all {eg. the Casimir pendulum: 
a finite plane inclined at an angle above an infinite surface) . The relation between the optical approach and the 
"semiclassical" approach is discussed further in Section HI. 

B. The sphere and the plane 

Next we analyze a problem with non-planar conductors — typical of real experimental configurations |f| — a sphere 
of radius R separated by a distance a from an infinite plane. In Ref. [g we tested the optical approximation by 
computing the Casimir force between a sphere and a plane up through the fourth reflection. We showed that the 
optical approximation is in very good agreement with the numerical results of Ref. l4!| for a/R ^ 1. In fact the 
numerical results presented in Ref. [g suffered from an insufficiently accurate numerical integration algorithm. The 
results presented here supercede Ref. |g| and show that the optical approximation is even more accurate than we 
originally claimed. For example, the optical approximation and the numerical data differ by only 30% at a/R « 5. 
Here we explain in detail how to compute the first and second reflection contributions. The relevant paths are shown 
along with some other aspects of the geometry in Fig. |2] For each reflection we must compute a) the optical path 
length, £r{x), b) the enlargement factor, Ar{x), and c) the domain of integration Vr for which r-reflections are 
possible. The Vr are subsets of the domain V above the plane and outside the sphere. 

1. Optical path lengths, £r{x) 

Finding the r-reflection optical path from x back to x, £r{x), is elementary in principle. One just draws straight 
lines from x to a given surface, from the arrival point on this surface to another surface, and so on, returning after 
r reflections to the original point x. One then moves the points of reflection on the surfaces until one reaches the 
minimum total length (an elastic string would do the job). The minimum length path suffers specular reflection upon 
each encounter with a surface. In all but the simplest geometries this problem must be solved numerically. However 
it is a problem amenable to extremely quick numerical solution: it is easily defined and the minimum is unique (at 
least for convex surfaces). This procedure also defines the points of refiection, Xr,i{x), Xr,2{x), etc. The first reflection 
paths from the sphere (Is) and the plane {Ip) and the two reflection path arc shown in Fig. 0. 
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FIG. 2: Geometry and reflections for a sphere and a plane. The regions and geometrical constructions are defined in the text. 

2. Enlargement factor, Ar (a;) 



The enlargement factor for the closed path beginning and ending at x is a special case of the general enlargement 
factor, Ar(a:,a:') for propagation from x to x\ well known in optics[3,|S|. In another guise, it is also well known to 
field theorists: Ar(a;, a;') is just the van Vleck determinant arising from the Gaussian fluctuations of the action about 
the classical r-reflection path from a; to a;'. In Section III, where we discuss the origins of the optical approximation, 
we show that the evaluation of the determinant gives the standard optics definition, 

A.(a:,x') = ^ (11-6) 

From this definition it is clear that in order to obtain Ar(a;, x') one must follow the spread of an infinitesmal pencil of 
rays of opening dfJj. from their origin at x, along this path, and measure the spread in area dA^i when they arrive at 
x' . Having aheady identified the points of reflection in lllB H it is relatively easy to compute A(a;, x'y\x,i=x numerically 
by tracing the paths of a few nearby raysQ- It is also possible to solve this problem analytically. Here we present the 
analytic solution for the first and second reflections from a sphere and plane. Beyond this level, it is probably more 
efficient to proceed numerically. One reflection from the plane is trivial: A(a;,a::) — l/£l{x) One reflection from the 
sphere is simplified by a) normal incidence, and b) x = x' . The second reflection can be simplified by regarding it as 
a single reflection from the sphere starting from x and ending at the image x of the original point x in the plane. In 
that case we need Ai(x,x). Consider the path from x to the sphere at the point Q, and then to x' . To obtain A one 
must follow the wavefront radii of curvature along the ray. We consider a ray that impacts the sphere at an angle 9 
to the normal. It travels a distance ui before and (T2 after, with (, — ai+ (T2. These variables are defined in Fig.j^^a). 
Consider a pencil of rays originating at x, spanning two infinitesmal arcs of angular widths rf0i,2 along perpendicular 
directions. Let dx'i and dx'2 be the associated arc lengths observed at x' . Then 
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Since both the initial ray and the sphere have equal radii of curvature we have the freedom to choose the directions 
defining d^i and d(j)2- We choose "latitude" and "longitude" as follows: Latitude is the direction perpendicular to the 
plane formed by x, x' and the center of the sphere (see Fig. |3fa)). Longitude is the direction in the plane. Consider 
the pencil of rays of varying latitude as shown in Fig. |2Ib). The variables are defined in the figure. It is easy to see 
that 



and considering that da = 2dfi cos 9 - 



dxi = aid(j)i + a2da, 
bi and Rdp — aid4>i, we find 
da 2(7i cos 9 
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FIG. 3: Geometry for reflection in a sphere, (a) The ray from x to x' reflecting at Q. Nearby rays originating at x and lying in 
the plane vary in longitude. Nearby rays out of the plane vary in latitude, (b) Variables for the calculations of the enlargment 
factor associated with latitude. The xx' plane has been projected along the vertical. A nearby ray originating at x heading out 
of the xx'Q plane by an angle (jii is shown. This ray reflects from the sphere at Q{4>\). The angle subtended by x and Q{<j)i) 
from the center of the sphere is /3. The angle formed by the vector from the center of the sphere to x and the ray from Q{<j)) 
to x'i{(j)i) is a. In the diagram the distance ci and the angles a and 01 are modified by factors of cos 5 due to the projection. 



and hence 
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The same calculation appUes for the longitudinal displacement except that the relation between d(/)i and dP is replaced 
by RdP = aid(f)2/ cos 9 and da = 2d(i + d0i, with the result 
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Putting these formulas together wc find for a single reflection on the sphere (the subscript s indicates reflection from 
the sphere) with angle of reflection 
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(11.12) 



Note that A(x,x') is symmetric with respect to the interchange of x and x' as it must, because the propagator 
possesses this symmetry. For the first reflection from the sphere we have cos^ = 1 and cti = (T2 = £/2, so 



Ai,(a;) = 



(£ + £y2R)^ 



(11.13) 



and, as mentioned above, the enlargement factor for the second reflection (on the sphere and then on the plane or 
vice versa), is given by the first reflection from x to its image x in the plane, A2(a;) = Ais(a;, x). A similar approach 
to higher reflections would require further analysis. The original wavefront leaving x is spherical. The first reflection 
from the sphere produces a new wavefront with, in general, two unequal radii of curvature. When next incident upon 
the sphere, the asymmetric wavefront will be transformed in a manner yet to be described. The ease with which 
Ar{x) can be computed numerically makes this unnecessary. 



3. Domain of the r^^ reflection, Vr 



The next step is the integration over the domains appropriate to each reficction. The first reflections give rise to 
cutoff dependent but a- independent contributions which must be analyzed at this point. Consider the first reficction 



from the plane. The appropriate domain is all of space except the interior of the sphere Vs and the region shadowed 
by the sphere "Dpi (see Fig.EJ- The integral can then be written as the difference between the integral over all space 
and the integral over Vg IJ "Dpi . The integral over the whole space is the divergent constant discussed in III Al which 
does not contribute to the force. It is to be ignored in the following. So the correct domain for the first reflection 
from the plane is the region, 'Ds[J'Dpi, in the shadow of the sphere, and the sign is to be reversed. Similarly, the 
integral of the first reflection on the sphere must be performed on the domain consisting of the whole space minus 
the interior of the sphere (Vs) and the region below the plate (T>p). The irrelevant divergence is given by the integral 
over all the space minus the interior of the sphere and the finite, a-dependent part, which contributes to the Casimir 
force, is given by the negative of the integral over the region, Vp, below the plane. So the correct domain for the first 
reflection from the sphere is Vp and the sign of the contribution is reversed. Hence we can write 

^^ f j3^ 1 ^c f 3 A}f (x) 



The second reflection gives a finite contribution to £. The path length, i2{x), never vanishes so there are no diver- 
gences at short distances, and the integrand, A2 {x)/£^{x), falls rapidly at large distances. The result is typically 
approximately 90% of the total result. Higher reflections can be analyzed in a similar fashion. The integration do- 
mains become progressively more restricted. For example, three reflection paths that reflect twice from the plane 
and once from the sphere do not exist in the shadow of the sphere (2?pi) nor in the darkly shaded regions in Fig. 13 
determined by the geometrical construction indicated by the dashed lines. It is not hard to carry out the constructions 
and calculations necessary to construct the optical approximation for the sphere and plane to any required order. 

4-. Discussion of numerical results 

In Ref. jy] we presented initial results on the optical approximation for the sphere and plane. Here we present final 
results (see Fig.QJ, discuss them in more detail, and compare them with the results of Ref. U and with the proximity 
force approximation (PFA). In presenting our results we display the sum of all the reflections (even and odd) up 
to (and including) the fourth. Since the energy must approach the parallel plate limit as a — > we can estimate 
the error in neglecting higher reflections in this limit. The error in neglecting the fifth and higher odd reflections 
is a -1-3.8% excess (because the sign of the odd reflections contribution is opposite to that of the total energy) as 
a ^ 0. Neglecting the even reflections (6th, 8th, etc.) as a ^ gives an error of —1.8%, negative because these 
contributions have the same sign of the total energy. Altogether the sum of the first four reflections overestimates 
the energy by 3.8% — 1.8% = 2% as a ^ 0. To illustrate this estimate of accuracy we have plotted our results as 
a band 2% in width in Fig. 2| Since the fractional contribution of higher reflections decreases with a, we believe 
this is a conservative estimate for larger a. Obviously, calculating the higher reflections will reduce this uncertainty 
interval, leaving eventually only the error due to diffraction which we are not able to estimate. The proximity force 
approximation has been the standard tool for estimating the effects of departure from planar geometry for Casimir 
effects for many years [l8|. In this approach one views the sphere and plate as a superposition of infinitesimal parallel 
plates: 

where d{x) is the distance from the plate to the sphere at a point x on the surface S. This formulation is ambiguous 
since the surface S could be taken to be either the sphere or the plate. Whichever surface is chosen, the distance is 
measured normal to that surface. The ambiguity is useful since it gives a measure of the uncertainty in the PFA. In 
either case the relevant integrals are easily performed. For the plate we obtain, 

cplatc TT^hcR 1 m1«^ 



while for the sphere we obtain 



^hcR 



a „ a / , a, 



In the limit a/ R —t both estimates agree to lowest order. The a ^ limit is usually called the proximity force 
theorem and has been much discussed over the years. It is usually stated as a result for the Casimir force in the limit 



a/ R -^ 0: j'^pFA ^ 2ttR£/A = —TT'^hcR/720a^ (where £/A is the Casimir energy per unit area for paraUcl plates). 
This hniit provides a convenient parametrization of the Casimir force when a is not so small, 



■'^R'720a3 



(11.18) 



Modern experiments are approaching accuracies where the deviations of f{a/R) from unity may be important. The 
accuracy of the PFA beyond the a/R — > limit is unknown, and the two different versions give different 0{a/R) 
corrections: 
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(11.19) 
(11.20) 



An important application of the optical approximation is to obtain a more accurate estimate of f{a/R). The optical 
approximation to the Casimir energy and the data of Ref. |j| both fall like l/a^ at large a. In fact both are roughly 
proportional to 1/a^ for all a. In contrast the PFA estimates of the energy falls like 1/a^ at large a and depart from 
the Gies et al data at relatively small a/R. For purposes of display we therefore scale the estimates of the energy by 
the factor —lAAOa^/Tr^RHc. The results are shown in Fig. 0] At large a/R the optical approximation has the same 



1.2 



1.0 

o «; 0.8 

"^ CO 



0.6 



0.4 



0.2 



















i 


Tr 


A A 






* 




■^ — *-= 


-*- 


— — 


^^-^^ A 










A 










S. 










\ 

X 
X 

V 


A- 


1 


1 




1 1 1 1 1 


• 



1/256 



1/64 



1/16 



a/R 



1/4 1/2 



FIG. 4: Casimir energy for a sphere of radius R a distance a above an infinite plane. 1440a S/ty Rhc is plotted versus a/R. 
The stars with error bars are the data of Ref. Q. The thick sohd curve is the optical approximation through the fourth 
reflection. The width of the curve indicates our estimate of the error in the optical approximation from neglect of the odd and 
even reflections with n > 5. The dashed curve is the plate-based proximity force approximation. The triangles are the results 
we published in Ref. Q, which are superceded by this work. 



scaling behavior as the data and differs from Ref. 4] by no more than 30% at the largest a/R. At small a/R, given 
our estimate of the accuracy of the optical approximation, we find that 
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(a/R) = 1 + O.Q5a/R + O {{a/Rf 



(11.21) 



which must be compared with the predictions of PFA eans. ljll.19p and ljll.20|) . In Fig. Owe show the contributions 
to the optical approximation of the different reflections we have computed. As expected the dominant contribution, 
always greater than 90%, comes from the second reflection. The fourth reflection contributes about 6% for a/R <C 1 
and less as a/R increases. The contributions of the first and third refiections are very small for all a/R. A relevant 
result, confirmed by the analytical analysis on the energy momentum tensor (within the optical approximation), 
the subject of the second paper of this series, is that the asymptotic behavior of 8 as a/R ^ 1 predicted by the 
optical approximation is oc 1/a^. This is in contrast with the Casimir-Polder law which predicts 1/a* at large a the 
discrepancy is to be attributed to diffraction effects. 



tJO 



1.0 



0.8 



Q^ 0.6 



I 0.4 



0.2 



2nd reflection 



1st and 3rd reflection 
from plate 



4nd reflection 



0.01 




1st and 3rd reflection 
from sphere 



FIG. 5: Contributions of specific reflections to the optical approximation. 



Casimir Pendulum 



In this section we treat a problem for which the exact answer is unknown. The configuration is shown in Fig. 1^1 The 
base plate is taken to be infinite. The upper plate is held at its midpoint a distance a above the base plate. The width 
of the upper plate is w and its depth, d (out of the page) , is assumed to be infinite. We define the Casimir energy 
per unit depth, e = £/d. 9 is the angle of inclination of the upper plate. It will be convenient to use z = ^wainO 
as a variable as well. It is also possible to view this configuration as a finite slice between £i = a/ sm9 — w/2 and 
£2 = a/ smO + w/2 of a wedge of opening angle 9. In this section we will discuss both the Casimir energy and the 
"Casimir torque" 1 '^ — hjmi psr unit depth. We are aware of two ad hoc approximate approaches to this problem. The 




FIG. 6; Geometry for the Casimir Pendulum. 

first is the PFA which treats each element of the system perpendicular to the lower plate as an infinitesmal parallel 
plate Casimir system. It is easy to show that 



£PFA 



TT^ hew COS 9 

1440 a3 



9 . 2 
w sm I 



4a^ 



n'^hca Vw^ 



Az' 



1440 (a2-z2)2 



which gives a torque. 



vpFA{a,w,z) 



n'^hc az{w'^ 
720 



32^) 



(a2 - ^2) 



2^3 



(11.22) 



(11.23) 



where the minus sign denotes that the torque is destabilizing: z = is a point of unstable equilibrium. As in the case 
of the sphere and the plane, the PFA is ambiguous. A more symmetric treatment of the two planes in the present 
geometry would integrate over the surface that bisects the wedge and take the distance normal to that surface. The 
result is the replacement of cos9 by cos'^{9/2) in eq. (|II.22|I and a similar modification of the torque. A second 
"approximate" treatment of the Casimir Pendulum can be extracted from the known exact solution for the Casimir 
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energy density for the "Dirichlet wedge" |l9j , which consists of two semi-infinite plates with opening angle 6 meeting 
at the origin. One can obtain an estimate of the energy for the pendulum by integrating the energy density over the 
two dimensional domain bounded (in polar coordinates, (p, i^)) by < (/) < and £i < p < £2- This approach takes no 
account of the modification of the energy density due to the finitenss of the upper plate. Furthermore it is inherently 
ambiguous because the energy density for a scalar field is only defined up to a total derivative. The calculation in 
Ref. !l9| used the conformally invariant stress tensor. One would obtain a different answer if one used, for example, 
the Noether stress tensor. In light of these difficulties, we do not pursue this approach further. To compute the 
optical approximation we need the enlargement factor, the lengths of optical paths, and the integration domain, Vr- 
Since all the conducting surfaces are planar, the enlargement factor is trivial in this case, Ar{x) —> l/i'^{x). The path 
lengths are also easy to compute. The only non-trivial step is the determination of the integration domains. As in the 
case of parallel plates, the odd reflections do not contribute to forces or torques for the Casimir pendulum. Instead 
they sum to a cutoff dependent constant associated with each plate. Any odd reflection path "turns around" with a 
reflection at normal incidence from one plate or the other. Consider all the points, x, which are the origins of paths 
that turn around at a given point P on either surface. These paths are shown, for the case where P is on the lower, 
infinite plane, in Fig. [3 They comprise one reflection paths lying on the interval PQi, three reflection paths lying on 
the interval Q1Q2, etc. The contributions to e from these intervals integrates to the same result as the integral over 
z for odd paths in the case of parallel plates. It is independent of a, w, and 9 and can be set aside. The fact that the 
enlargement factor is trivial is crucial for this argument. 




FIG. 7: Odd reflection paths for the Casimir Pendulum. 



1. Even optical path lengths, £r{x) 

The analysis of paths that reflect an even number of times makes use of simple geometrical concepts. For any point 
X = (p, 4>) , we define an infinite sequence of images in the upper and lower planes as shown in Fig. |S1 ignoring for the 
moment that the upper plane is finite. The images below the lower plane are denoted xi,X2, ■ ■ ■ and those above the 
upper plane are denoted a;i,a;2, .... In sequence, xi = Rx, xi — Rx, X2 — Rxi, X2 — Rxi, etc. with R denoting 
reflection in the lower plane and R denoting reflection in the upper plane. All of the images lie on the circle of radius 
p about the origin. The length of the 2n reflection path can easily be seen to be given by 

i2n{x) = \\xn - x„ 1 1 = V sin^ ^6 (11.24) 

independent of 0. Substituting into eq. (|I.1|I we obtain the expression for the Casimir energy per unit depth. 



-opt 



he -^^ 

167r2 ^ 

n— 1 



sin'' 126 



1 

dp^ e(p2n) 



(11.25) 



The step function Q(T>2n) vanishes when the point x is not in the domain where 2n-reflection paths are possible. As 
we show in the following section, for n large enough, P2n -^ 0- -/V,„ax, the upper limit on the n-sum, is the largest 
value of n for which any 2n-reflection paths exist. 



2. Domain of the 2n^^ reflection, 'V2n 



The domain in which the 2n-reflection path exists is determined by the constraint that the points of reflection at 
the upper plate must lie between £1 and £2, the inner and outer radii that deflne its boundaries. Note, of course, that 
(.2 > ii- Although the calculation is elementary, it is tricky, so we only quote the results. The constraints depend on 
whether n is even or odd, so we summarize them independently. 
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(p,-^) 



FIG. 8: Images of the point x in the Casimir Pendulum configuration. The dashed hues have the same length as the r — 2,4, 6, 
reflection paths. 



• 71-odd When n is odd, the integration domain I?2n is defined by the inequalities: 

£i cos (j) < p cos n9 < £2 cos((n — 1)9 + (j>) 



(11.26) 



where the lower limit ensures that the innermost reflection occurs at p > £1 and the upper limit ensures that 
the outermost reflection occurs at p < ^2- The inequality cannot be satisifed for any p or 0(< 9) unless 



£2Cos{n-l)9> £1 

When eq. IJII.27P is satisfied the contribution of the 2n*'' reflection (n-odd) is 

he cos^ n9 f |t tan6' — 1^ (tann^ — tan(n — 1)0) for £2 cos ri6' > ^1 cos 6 

for £2 cos nO < £\ cos ( 



.odd 
-opt n 



327r2 sin-* n9 



{i2cos{n-i)e-eif 



£j£| cos(n-l)e sin(n-l)6 



• ?i-even When n is even, the integration domain I?2n is defined by the inequalities: 

£1 cos(6' - (p) < pcosn9 < £2 cos{{n ~ 1)9 + (f) 



(11.27) 



(11.28) 



(11.29) 



where, as before, the lower limit ensures that the innermost reflection occurs aX p > £1 and the upper limit 
ensures that the outermost reflection occurs at p < ^2- The inequality cannot be satisifed at all unless 



I2 cos(n -1)9 > £1 cos 9 
When eq. Ijll.30|) is satisfied the contribution of the 2ii}^ reflection (n-even) is 



-opt n 



he cos^ n9 
327r2 sin'^ n9 



■ tan ( 



■h (tann^ — tan(n — 1)6*) for £2 cos n9 > 



(£2 cos(7i-i)e-f 1 cos e)^ 

I^£^ cos(n — l)9 s'lnnO cos ^ 



for £2 cos n9 < £\ 



(11.30) 



(11.31) 



The torque is obtained by differentiating with respect to 9 at fixed a and w, remembering that £\ and £2 depend 
GVi 9. Of course the 6'-derivative of eqs. (|II.3ip and (|II.28I) are complicated and need not be written down explicitly. 
The resulting expressions for v'^t^ ^ ^^^d ^opt"n must be summed over n subject to the constraints in eqs. I|II.28|) and 
(|11.31|l . This sum must be performed numerically. The results are discussed in the following subsection. 



Discussion 



Fig.|51shows the pendulum Casimir energy as a function of z for a ~ 1 and several values of w. The weak dependence 
on z at small z is to be expected. So is the divergence as 2; ^ 1 which we do not show in the figure. When the plates 
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FIG. 9: Casimir energy and torque in scaled units for a Casimir pendulum of width w = 1.5a, 2.5a, 5a and 10a. Positive values 
of the scaled torque are destabilizing. 

touch at z = 1, the Casimir energy of perfectly sharp, perfectly conducting plates would in fact diverge, as would the 
Casimir torque. The optical approximation for the pendulum turns out to be very close to the plate based PFA. It is 
convenient therefore to scale out a factor of — 1440a^/?ic7r^w (see eq. (jll.22|) ') when displaying or results for the energy 
per unit depth, Copt, and a factor — 720a^/?ic7r^w^ when displaying results for the torque per unit depth. The energy 
and torque are plotted for representative values oiw (w = 1.5a, 2.5a, 5a, and 10a) as a function of zja in Fig. El The 
plots are shown only up to zja — 0.5. Above zja ~ 0.5 both grow rapidly and diverge at z — a. The PFA gives an 
excellent approximation for the pendulum over a wide range of the parameter space. This can be seen by examining 
the ratio of eopt/epFA as shown in Fig. ^| The reason behind this success is that the second (optical) reflection is 
proportional to the PFA result (eq. III. 2211 ) for all z. The constant of proportionality is the familiar 90/tt^ — 0.924. 
The sum of the higher even reflections combines with the second to equal the PFA at z = and drops away slowly with 
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FIG. 10; The ratio of the optical approximation to the PFA. (a) for a pendulum of different widths as a function of z/a. The 
breaks in the curves for w = 1.5a and 2.5a occur when only the second reflection can contribute. The w = 1.5a curve ends at 
the z — 0.75a when z = w/2. (b) The contributions of the 4th and 6th reflections for w — 2.5a. The 2nd reflection contributes 
0.924. . . independent of z. This is the only case we found in which the optical approximation gives an energy smaller than the 
PFA. 



increasing z/a. So the optical estimate coincides with the PFA at z = and drops slowly with z. The contributions 
of the first few reflections are shown for w — 2.5a in Fig. ^| A careful study of Fig. |51 reveals one peculiar feature of 
the optical approximation which is interesting in principle, if not in practice. The torque does not vanish at z ~ 0, 
and therefore the Casimir energy has a cusp at z = 0. The section of the graph near z = is enlarged in Fig. 1111 to 
make the effect clearer. The torque starts negative (stabilizing) at very small z before it changes sign and becomes 
destabilizing as z increases. This is an effect of the finite extent of the upper plate. It vanishes like a/w at large 
w. We are not able to determine whether this is an artifact of the optical approximation's neglect of diffraction, or 
whether it is a real effect at small z. We know of no theorem that precludes a cusp in the energy at z = 0. We will 
have to wait until the effects of diffraction can be included before analyzing this further. 
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FIG. 11: Torque in scaled units for a Casimir pendulum of width w = 1.5a, 2.5a, 5a and 10a, for z ~ 0. Positive values of the 
scaled torque are destabilizing. 

III. ORIGINS OF THE OPTICAL APPROXIMATION 

Most studies of Casimir energies do not consider approximations. Instead they focus on ways to regulate and 
compute the sum over modes, ^ ^fi(i;|l8|. These methods have proved very difhcult to apply to geometries other than 
parallel plates. The main reason for this impasse lies in the requirement of an analytic knowledge of the spectrum. 
Finding the spectrum of the Laplace operator for non-separable problems is not merely a technical difhculty, it is 
more one of principle. In fact there are strict relations between this problem and those of chaotic billiards theory. 
The existence of an exact solution for the Casimir problem with a non-trivial geometry would imply the existence of 
an exact solution for the same family of quantum billiards and hence of classical billiards. 30 years of work on the 
ergodicity of classical billiards and the implications for the density of states in the corresponding quantum billiards 
suggest this task is hopeless (see jlSj). On the other hand, one could argue that a numerical solution of the Dirichlet 
problem could easily give the spectrum to some high, but finite, accuracy, and this could be used to compute the 
Casimir energy. However the Casimir energy diverges - the leading non-trivial divergence in N dimensions is of order 
A^. The Casimir force between rigid bodies is known to be finite, and one could hope to compute it by introducing 
a cutoff, computing the energy at nearby separations, a and a + da, taking the difference, £{a + da, A) — £{a. A), and 
finally taking A — > oo. However such numerical problems are hopelessly unstable. The force, indeed, is given by the 
small oscillatory ripple in the density of state numerically shadowed by the 'bulk' contributions which give rise to 
distance-independent divergencies. So we focused our attention on ways to get approximate solutions of the Laplace- 
Dirichlet problem which are apt to capture the oscillatory contributions in the density of states, providing physical 
insights and accurate numerical estimates. We have not found any previous use of ideas from classical optics. In this 
section we give a derivation of the optical approximation based on a path integral representation of the Helmholtz 
Greens function. Schaden and Spruch have developed an approximation for Casimir energies |ll| using Gutzwiller's 
semiclassical treatmentfl^ of the density of states. It is misleading to call the approach of Ref. [Ill "semiclassical" 
because, as can be seen for example from eq. ()III.1|I . the only h in the Casimir problem for a massless field is the 
multiplicative factor in ^huj. However, since the authors of Ref. [llj use the term following Gutzwiller, we will continue 
to refer to their approach as "semiclassical" . This work differs in important ways from ours and in general is not as 
accurate, however the relationship between the two approaches is interesting, and is explored later in this section. 

A. Derivation 

We begin with the well-known definition of the Casimir energy in terms of a space and wavenumber dependent 
density of states|20|, p^x, k), 



£vm 



dk d x—h uj{k)p{x,k), 
Jv ^ 



where ijj{k) — ci/fc^ + y?, and the density of states p{x, k) is related to the propagator G{x' , x, k) by 



2fc 
p{x, k) — — Ini G{x, X, k). 

n 



(IILl) 



(IH.2) 
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Since we are considering a scalar field, G is the Greens function for the Helmholtz equation. We choose G to be 
analytic in the upper-half fc^-plane (or equivalently take k^ to have a small positive imaginary part). The tildes on 
p{x,k) and G{x,x' ,k) denote the subtraction of the contribution of the free propagator, Gq{x' ,x,k). The Casimir 
energy depends on the boundary conditions obeyed by the field -0 and on the arrangement of the boundaries, S = dV 
(not necessarily finite), of the domain V. From the outset we recognize that £ must be regulated, and will in general 
be cutoff dependent, as discussed in the Introduction. We will not denote the cutoff dependence explicitly except 
when necessary, p is the familiar density of states associated with the problem 



so that G satisfies the equation 
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(III.3) 



(III.4) 
G{x', X, k) = G{x', X, k) - Go(a;', x, k), (III.5) 



where Gq is the free scalar propagator in the absence of boundaries. The spectral representation expresses G as a 
sum over a complete set of eigenfunctions tpn with eigenvalues fc„ 

G(x,x,/c) = _^^^— ^^— -. (III.6) 



Notice that since the problem (|III.3|I is real we have chosen a complete set of real eigenfunctions and removed the 
usual complex conjugation from (|III.6|I . We can regard this problem as the study of a quantum mechanical free 
particle with h ^ \, mass m — 1/2, and energy E = k'^, living in the domain T> with Dirichlet boundary conditions 
on dT). Dirichlet boundary conditions are an idealization of interactions which prevent the quantum particle from 
penetrating beyond the surfaces S. This idealization is adequate for low energies but fails for the divergent, i.e. 
cutoff dependent, contributions to the Casimir energv|21|. As we have already seen in Section II, the divergences 
can be simply disposed of in the optical approach, and the physically measurable contributions to Casimir effects are 
dominated by fc ~ 1/a, where a, a typical plate separation, will satisfy 1/a <C A where A is the momentum cutoff 
characterizing the material. So the boundary condition idealization is quite adequate for our purposes. Following this 
quantum mechanics analogy we introduce a fictitious time, i, and consider the functional integral representation of 
the propagator p^. The space-time propagator is 

/°° (7P 
G{x',x,^)e-'^\ (111.7) 

where E = fc^. Since G is analytic in the upper half fc^-plane, it is evident that G{x',x,t) = when t < 0. The 
inverse Fourier transform reads 

/•oo 

G(x', x,k)=i dt e*'='*G(a;', x, t). (III.8) 

G{x',x,t) obeys the free Schrodinger equation in D bounded by S. It can be written as a functional integral over 
paths from x' to x with action S{x' , x,t) ^ j J dtx^ . The optical approximation is obtained by taking the stationary 
phase approximation of the propagator G in the fictitious time domain. Hence we assume that the functional integral 
is dominated by the contribution of classical paths between x' and x. These are straight line paths, reflecting r times 
from the boundaries, and traversed at constant speed, v — ir{x' , x)/t. where ir{x' , x) is the length of the path. Then 
the optical approximation to the propagator is given by, 

Gopt(x',x,i) = ^i?,(x',x,t)e'^'-(^'^'*). (III.9) 

r 

The action is 

5.(x',x,t) = ^li^ (III.IO) 
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/ Q2(2 \l/2 
Dr{x\x,t)^dci[^—r-] {lll.ll) 



and D is the van Vleck determinant 

\^x^^Xj 

This approximation is exact to the extent one can assume the classical action of the path Sr to be quadratic in 
x',x. This is the case for flat and infinite plates. Thus the non-quadratic part of the classical action comes from the 
curvature or the finite extent of the boundaries, which we parameterize generically by R, d'^S/dx^ ^ 1/Rt. Hence, 
in a stationary phase approximation Sx ^ ^/i and the corrections are of order 6^S ^ O (Vi/i?). Back in fc-space the 
corrections hence will be O (l/kR), and the important values of k for the Casimir energy are of order 1/a, where a 
is a measure of the separation between the surfaces. Thus the figure of merit for the optical approximation is a/ R. 
At the moment there is no good way to estimate the order in a/ R of the corrections to the optical approximation 
(possibly fractional, plus exponentially small terms). Certainly some of the curvature effects are captured by the van 
Vleck determinant, and as we saw in section FlI Bl for the sphere-plare problem, the optical approximation works in 



practice out to a/R ~ 1. This is topic for further investigation. Eq. (|III.9|I is, in fact, the usual approximation of ray 
optics, the van Vleck determinant being precisely the enlargement factor of classical optics, as we now show. Since 
dir{x' , x)/dx' = n' and d£r{x' , x)/dx = — n, where n and n' are the unit tangent vectors to the path in the points x 
and x' , we have 

Dr{x',x,t) <xdetln,n'^+£rj-j-] , (III.12) 

We perform the analysis in three dimensions. Other values of N are analogous, and we quote the general result at 
the end. The matrix -a-^ is 

^ti®t; + ^t2 0t'„ (III.13) 

where ti.2 and t'j^ 2 ^^^ orthonormal tangent vectors perpendicular to n and n' respectively and with them form two 
orthonormal bases centered in x and a;', and d(j)i/dx[ is the derivative of the angle subtended at the point x when we 
shift the point x' along the direction t[. Taking the determinant is now easy: it is the product of the three eigenvalues 
of the matrix, but given the fact that {n, ti,t2} (and their primed correspondents) are an orthonormal triple these 
are just {l,£rd(j)i/dx[,ird(l)2/dx2}, so 



Dr{x',X,t)^[ilf^^ 



1/2 

(III.14) 



The coefficient of proportionality is independent of the path^ r and must depend on t in such a way that for the direct 
path we obtain the free propagator. Therefore, 



(-l)"^ f,N-idn 



where we have returned to A^-dimensions. We have introduced the factor (— 1)*^ to implement a Dirichlet boundary 
condition. In the case of a Neumann boundary condition, this factor would not be present. Although we did not label 
dft/dA with an index r, it should be clear from the derivation that it does depend on the path r. Putting all together 
we find the space-time form of the optical propagator to be 

ftp.«M)^i:^5yife(c-'^) "%•'='« (....«) 

When dealing with infinite, parallel, flat plates this approximation becomes exact. For a single infinite plate, for 
example, the length-squared of the only two paths going from x to x' are 

(III.17) 
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^ We are not discussing the Maslov indexes other than the (— 1)*^ here. If the ray r would touch a caustic it would be necessary to 
introduce the appropriate phase factor. 
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where x is the image of x. Both are quadratic functions of the points x, x' and the optical approximation is indeed 
exact. In order to calculate the density of states we must return to fc-space. G{x',x,k) is obtained by Fourier 
transformation (see eq. Hill. 8(1 ). and can be expressed in terms of Hankel functions, giving us the final form for our 
approximation 

= ^Grix',x,k), (III.18) 

r 

where N is the number of spatial dimensions, A^ is the enlargement factor 

Ar{x',x) = ^ (111.19) 

and we have suppressed the arguments x and x' on £,. and A^ in l|III.18|) . This can be thought of as a particular case 
of the general results in Ref. J23| • 

For A^ = 1 and iV = 3 the Hankel function reduces to an exponential. For example, when A^ = 3 we find 



A 



1/2 



Gr{x',x, k) = (_i)"^^!_e*'=^'-. (III.20) 

47r 

However, had we attempted a stationary phase approximation directly in /c-space we would have obtained an expo- 
nential for any N, 

Gsomicl{x',X,t)= Y. Dr{x',X,k)e''''^^^'^\ 
paths r 

which does not reduce to the exact expression in the limit in which we have only infinite, non-intersecting (hence 
parallel), flat planes, because in k space it is not a gaussian problem even for quadratic t"^. This is an important 
advantage of applying the stationary phase approximation in the time domain where it leads to the optical approx- 
imation. Also we believe the optical approximation to be a more favorable starting point for considering systematic 
corrections to the stationary phase approximation uniformly^ in 1/R. The expressions pil.9|l and l(III.18|l are the 
first term of a systematic expansion of the propagator in 1/kR. For gently curved geometries we expect them to 
provide a good approximation, the final test, in absence of exact solutions, coming only from comparison with the 
experiments. The corrections come from two different (but related) effects |25j : a) we have to expand the function 
Sr{x',x) in the exponential to include cubic (and higher order) terms and b) we have to include other stationary 
paths of non-classical origin, like paths running all around the bodies one or more times (these can be considered as 
a non-perturbative, exponentially small correction to the propagator). Both phenomena are due to the curvature of 
the boundary surfaces and we go back to the previous estimate that the parameter controlling the accuracy of the 
our approximation is indeed 1/kR (wedges and discontinuities must be considered as regions in which R —f and 
the expansion is somewhat different) . Two intertwined branches of wave optics have dealt with finding corrections to 
the geometric optics predictions for curved boundaries. The first 8,'l5l deals both with perturbative a) and nonper- 
turbative b) corrections to next to leading order in 1/kR of particular importance in the shadow region. The second 
deals with edges and holes in locally flat surfaces, originated by Sommerfeld's work J2g (see also [23 and references 
therein). Both must be considered relevant to future studies of Casimir forces, since high-curvature and flnite-size 
effects will soon be relevant in the next generation of precision experiments jla 1231 • Another phenomenon to be taken 
in account, even in the case of gently curved surfaces, the optical approximations fails when either x or x' are in the 
shadow region or we are in presence of a caustic, the set of points where the Hessian d^Sr/dxdx' has one or more 
vanishing eigenvalues |33- In these regions of the parameters (x, x') the gaussian approximation fails and one cannot 
ignore cubic terms in the action. There are various ways of treating this phenomenon, whose importance in wave 
optics|2J, [SJ as well as quantum mechanics p3l |25| is today clear. The most interesting prediction related to the 
presence of caustics (for what concerns us here) is the fact that a ray crossing a caustic acquires a non-trivial phase 
shift. This could possibly result in a change of the sign of the Casimir force for concave geometry. Unfortunately, the 
existing formalism does not seem to be easily translated into our language and more work is needed in this direction. 



The technique of passing to the Fourier transform to obtain uniform approximations is certainly not new in wave optics |24| . 
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The famous MRE of Balian and Bloch 32] is also intimately related to the optical approximation developed here. It 
is relatively easy to see that our approximation arises as the first term in a uniform 1/kR expansion for the propagator. 
Most of the effort in applying the MRE to Casimir energies has focused on the divergent terms associated with general 
geometrical properties of the bodies 133] or on the Casimir force at large distances where only the lowest reflections 
contribute. To our knowledge no one has been able to develop a useful expansion beyond the optical limit from the 
MRE. 

B. The optical Casimir energy 

The substitution of l|III.18(l into (|III.2p and then in (|III.ip gives rise to a series expansion of the Casimir energy 
associated with classical closed (but not necessarily periodic) paths 

paths r 

where each term of this series will be in the form of 

1 foo ojf r 

£r = -hlm dkuj{k)— d^x Gr{x,x,k). (III.22) 

Here the integration has been restricted to the domain 7)^ where the given classical path r exists. At this point it 
is useful to separate potentially divergent contributions from those which are finite. Because G is analytic in the 
upper half fc-plane, the k integration can be taken along a contour with Imfc > 0. The Hankel function i?„ (fc£) falls 
exponentially in the upper half plane, so the x integral converges absolutely and uniformly at fixed-fc unless there are 
X- values where lr-{x) can vanish. One can easily convince oneself that for smooth surfaces^ the only paths that reflect 
once on any surface can give vanishing path length. So for the moment, we put aside the first reflection and consider 
the cutoff independent contributions from r > 1. In that case we can interchange the k and volume integrals. The 
resulting fc-integral is also uniformly convergent. 

£. ^ ^ (1^ Re j^^ d-x^^l^^^ j^ dfc.(fc)^(2fc)-/-^<,_,(M.), (III.23) 

for r > 1. The /c-integral can be performed in general, but is particularly simple for the massless case, uj{k) = c/c, 

£r = ^C-^E^r (^) I d-X-^, (III.24) 



2^JV/2+l/2 \^ 2 J J-p^ £^^+3^/2' 

which is the Casimir energy associated to the optical path r > 1, and generalizes our fundamental result, eq. (|I.ip to 
dimensions other than three. The generalization to the massive case for A^ = 3 is given by 



vr+l^C^ [ ,Z^^r 



1/2 



Er = {-ir+'-^ d'x^—K^itltr), (III.25) 

which reduces to the N — 3 case of eq. (jIII.24p as /i -^ 0. It is worth nothing that for /i > the paths with length 
£ ^ 1//Li are exponentially damped. Now we return to analyze the potentially divergent first reflection. For simplicity 
of notation we specialize to A^ = 3 although the analysis is completely general. Let the boundary of D be the surfaces 
of a set of rigid bodies Bi, B2, ■■■, Bn- The divergent contributions come from the paths IBi that reflect once on any 
of the bodies Bi. To regulate possible divergences we insert a simple exponential cutoff in k. It is easy to see that 
our results are independent of the form of the cutoff. Then for a massless field, reflecting from body B, 

£iB = {-l)—:^ d^xAY^{x,x) dk e-''/^k^ sm{keiB{x,x)). (III.26) 

47r^ Jv^B h 



It suffices that the vector n normal to the surface is continuous, i.e. no wedges are present. 
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The fc-integration can be performed, 

Notice that for ^isA ^ 1 we reobtain the standard resuh, eq. I|I.1|I as we should. When £isA < 1 however the 
structure of the function changes completely. In particular the sign changes at £ibA = ^/i. There is a non trivial 
consequence of this fact: from eq. I|I.1|) one expects a positve divergence (r = 1 here) as ^ ^ 0, however the small (. 
divergence in eq. (|III.27|I is negative. This effect, that the cutoff dependent contribution to the Casimir energy density 
changes sign near the bounding surface, is well known and has figured centrally in recent discussions of Casimir energy 
densities 'id]. Of course the bulk contribution to the vacuum fluctuation energy comes from the zero-reflection term, 
which is positive. The negative surface correction is well known and has many physical consequences. For example it 
contributes to the surface tension of heavy nuclei 34] . 

To analyze the divergent first reflection, eq. IJIII.26P further, we need an expression for A(a;,a;) near a generally 
curved surface. This entails a small change in Ai^ (see eq. (jII.12p ') to take in account two different principal radii of 
curvature, say Ra and i?b (here x' = x so 9 = Q and (Ji = (T2 = ^/2), 

Substituting back into eq. pil.27|l and replacing (fx = dS{i)d£/2, where dS{i) = {e/2Ra + l)i£/2Rb + l)dS, and dS 
is the surface area element on the body, we get (up to finite terms arising from upper bounds on the integration in 
di) 

s^-^ldsj^^ diVii + e/2Ra){i + e/2Rt) ^^^^l~J[^J^7 ■ ("1-29) 

where we have suppressed the subscript IB. The ^-integration may be evaluated at large A to obtain an asymptotic 
expansion of the cutoff dependent terms in the first reflection, 

Eq. (|III.30p summarizes the cutoff dependent contributions to the Casimir energy in the optical approximation. As 
discussed in Section II, these terms do not contribute to the forces between rigid objects. Also they are trivial to 
isolate and discard from the calculation of forces. The form of eq. (|111.30|l invites comparison with the work of Balian 
and Bloch,^ on the asymptotic expansion of the density of states based on their Multiple Reflection Expansion. 
The MRE propagator includes not only specular paths, but also contributions from diffraction which also yield cutoff 
dependent contributions to the Casimir Energy. Scaling arguments indicate that terms up to at least the third 
"reflection" in the MRE are cutoff dependent. These higher divergences are omitted from the optical approximation, 
which is convenient since they do not contribute to Casimir forces in any case. The first few terms in the MRE 
expansion of the density of states are given by, 

so the leading cutoff dependent terms in the Casimir energy are* 

1 f°° 

£ ^ -he dkk pMRE(fc) 6-'=/^ 

^ Jo 
- -^^^A^-2i.^^A^/^^(^ + ]^)+^(A) (ni-32) 



The sign of the second term here is opposite that of Ref. j3^ because we are dealing with convex rather than concave geometries 
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Comparing with the optical result, eq. Ijlll.30|) we see that the first terms agree and the second terms differ by a 
factor of 3/4. Apparently our optical approximation to the propagator, despite its simplicity, captures the leading 
divergence and the order of magnitude of the subleading divergence^. 

C. Connections Avith other semiclassical approximations. 

Stationary phase approximations are not new in the study of Casimir energy both at zero and non-zero temperature 
[ill ISSf . These works certainly share with ours the attempt to switch the attention toward general properties and 
approximations to the Helmholtz equation. On the other hand, relying more or less heavily on Gutzwiller's trace 
formula, they suffer from two significant problems. First, they treat symmetric and nearly symmetric geometries in 
radically different ways, and fail to provide a natural deformation away from the symmetric limit (not to mention 
that they give the exact result for parallel plates only for odd number of space dimensions). Second, they require 
a certain amount of strongly geometry-dependent work (to calculate monodromy matrices for example) . We discuss 
both these problems further below. In order to study these points we will rewrite a given contribution £r (specializing 
to A'^ = 3 dimensions and suppressing the index r) as 

/•oo />oo 

£ = Im(-l)" / dk hck^e-^'^ / M J(^)e''=^ (III.33) 

where 

J{i)= [ d^xS{i-er{x,x))^^—^^^^^. (III.34) 

Our strategy has been to dominate the functional integral over paths from x back to x by the classical paths, then 
perform the k integral analytically, and to leave the integration over x for numerical evaluation. The standard 
"semiclassical" approach [ll|, [lal is to perform all spatial integrations by stationary phase including the one over the 
argument of the Greens function itself. This leaves a function only of k which can be integrated analytically. The fact 
that we can do the x integral numerically allows us to capture much more detailed information about the system. We 
will show this in detail in the following. To underline the differences, let us repeat briefiy the line of reasoning leading 
to Gutzwiller's trace formula. We start by writing an asymptotic expansion for k£ ^ 1. The asymptotic contributions 
to the ^-integral come 36| both from a) boundaries at im,iM (minimum and maximum length achieved by the path 
r) that is integration by part terms and b) integrable divergences in the function J{£) that is stationary phase (SP) 
points at £j. So that 

J{£) - J2 (Ai(fc)e*""' - Bn{k)e'^^'-) + ^ Q(fc)e^"^ (III.35) 

n>0 SP points j 

A,B,C are polynomial in fc, 1/fc and £jn,£M,(ij respectively. The Schaden and Spruch 11 1| approach based on 
Gutzwiller trace formula [l6| consists in taking only the stationary phase contributions, b), to the energy, the coeffi- 
cients Cj(fc)'s then being related to the 'monodromy matrix'^. These terms correspond to closed classical paths, for 
which the final momentum is equal to the initial one (the action is 5* ex ^). The stationary phase approximation re- 
quires the periodic orbits to be well-separated in units of wavelength. However as one approaches a situation in which 
one exact symmetry exists, the space, M^ in 3 dimensions, perpendicular to the closed orbit at a given point breaks 
into the product of two subspaces A(S) B, and £ is constant with respect to the B coordinates b. In the symmetric 
situation the SP points then form lines (or planes if more than one symmetry is present) parameterized by b. The 
problem can again be solved easily just by writing d^x oc dadb and factoring out the integral over (i& [llll37j leaving the 
integral over da to be evaluated by stationary phase approximation again^. However, when the symmetry is slightly 



One might think to claim more than order of magnitude success here. However it should be noted that for Neumann boundary conditions 
both terms in lIII.30t change signs while only the surface terms in 1III.321 changes sign. This is due to the fact that 2 'reflections' in the 
MRE expansion contribute to the curvature divergence as well and their sign is the same for Dirichlet or Neumann boundary conditions. 
To be precise the stationary phase integral is done on the directions transverse to the periodic orbit. The integration over the direction 
parallel to the orbit is eventually performed by means of a trick 1^ . 

The simplest example is that of a cylinder facing a plane. Then the periodic orbits are lines perpendicular both to the cylinder and 
the plane, b is parallel to the axis of the cylinder and a is the direction perpendicular to this. In the case of parallel plates both the 
directions a and b are symmetry directions so they both factor out and no stationary phase approximation is performed. In this case 
the former analysis gives an exact result, as is well known. 
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broken the length i acquires a smah h dependence and the integral over db can no longer be factored out. Moreover a 
naive stationary phase approximation in both dadb is not reliable because arbitrarily close to the breaking point, the 
dependence of £ on 6 is small and the Hessian matrix d'^(.(x,x) jdx^ has one (or more) very small eigenvalues in the 
old h directions. There exists |2a| a theory for Gutzwiller trace formula for approximate symmetries. However, we 
found that it is not easy to implement in the study of the Casimir energy for arbitrary surfaces. In the cases in which 



1.01 




FIG. 12: Comparison between the optical approximation (upper curve) and the "semiclassical" approximation of Schaden and 
Spruch (lower curve) for the sphere and plane. The scaled Casimir energy is plotted versus a/R. For the optical approximation, 
the sum of the first four reflections has been rescaled to go to unity as a — > 0. It is possible to show that in the limit a/R ^ 1 
the optical approximation and Schaden and Spruch's formula agree (in the flgure they both tend to OO/tt^ = 0.92...). The most 
notable and relevant discrepancies are in the derivative at small a/R. 

these problems can be avoided, like a sphere of fixed radius in front of a plane, and for a/R <C 1 the semiclassical 
approximation a la Schaden and Spruch provides quite a good approximation (see figure) since in the expansion 
(jIII.35|l the stationary phase approximation gives a much larger contribution than the integration by parts terms. 
A much stronger disagreement has to be expected if the sphere gets substituted by a plate of width w bent with a 
curvature of order R ^ w. Indeed the method of Ref. |llj| differs dramatically from the optical approximation for 
the case of a hyperboloidQ- This is not a diffraction effect but rather a 'precocious' breakdown of the semiclassical 
approximation which is cured by a uniform approximation of the kind we have described. 



IV. CONCLUSIONS 



We have proposed a new method for calculating approximately Casimir energies between conductors in generic 
geometries. We use a stationary phase approximation imported from studies of wave optics that we have therefore 
named the "optical approximation" . In this paper, the first of the series, we have outlined the derivation and applied 
it to three examples: the canonical example of parallel plates; the experimentally relevant situation of a sphere facing 
a plane; and the "Casimir pendulum", i.e. a conducting plate free to oscillate above an infinite plate, where the 
calculations can be performed analytically. In all of the above examples (except for parallel plates, where our result 
coincides with Casimir solution) the agreement with the Proximity Force Approximation is only to the leading order 
in the small distances expansion. The first order correction is found to be different. This is of particular importance 
in the example of the sphere and the plane because the first order correction in a/R (a is the distance sphere-plate and 
R is the radius of the sphere) will soon be measured by new precision experiments |29l | . The optical approximation 
turns the Casimir sum over modes into a sum over topologically different paths, and from this point of view can be 
compared with the Poisson summation formula, which has proved useful to derive semiclassical uniform expansions 
for very diverse problems 23, 39]. In the case of the Casimir energy, replacing the usual highly divergent sum over 
modes by a sum over topologically distinct optical paths has two, very significant advantages: first, we have been able 
to show that the divergences in the Casimir energy are contained in contributions of very simple, one-reflection paths 
and can be easily and unequivocally regulated and discarded; and second, the convergence of the sum over paths is 
very rapid. Instead of requiring an infinite number of eigenvalues with exquisite precision one needs but a few path 
contributions, calculated with little numerical effort, to give a very good approximation to the Casimir energy for 
important geometries. "Semiclassical" methods have been used previously in the study of Casimir effects and the 
connection between closed orbits and the finite part of the Casimir energy has been pointed out many times |l7l |. 
Our analysis shares with those the idea of shifting attention to approximations and to properties of the Helmholtz 
propagator. We have shown however that in order to obtain a correct low curvature approximation one has to use 
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a uniform approximation of the kind we proposed here. There is plenty of room to improve the approximation 
presented here, especially when the connection with Balian and Bloch's multiple reflection expansion is made explicit. 
In particular it is intriguing that the Casimir energy for the sphere-plane problem is a well-defined problem in a 
single variable, namely x — a/R whose limiting values for a; ^ 1 and a; S> 1 are famous jS."! l4(l| . One can hope that 
an analytic solution or an approximation good for the entire range z-values should be relatively easy to find. On 
the contrary it is an incredibly difficult problem and nobody has succeeded in finding such an exact solution or a 
valid approximation. In the next paper we will show how the same approximation for the propagator can give useful 
expansions for local operators like the energy-momentum tensor, which allows us to calculate the pressure the energy 
density and other properties of the constrained field. We will show what our approach has to say about the thermal 
corrections, which are a controversial matter in Casimir physics. In another paper of the series we will also perform 
the same analysis for a field of spin 1/2 and and for the electromagnetic field. 
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